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A realistic quantum many-body system, characterized by a generic microscopic HamUtonian, is accessible 
only through approximation methods. The mean field theories, as the simplest practices of approximation 
methods, commonly serve as a powerful tool, but unfortunately often violate the symmetry of the 
Hamiltonian. The conventional BCS theory, as an excellent mean field approach, violates the particle 
number conservation and completely erases quantumness characterized by concurrence and quantum 
discord between different modes. We restore the symmetry by using the projected BCS theory and the exact 
numerical solution and find that the lost quantumness is synchronously reestablished. We show that while 
entanglement remains unchanged with the particle numbers, quantum discord behaves as an extensive 
quantity with respect to the system size. Surprisingly, discord is hardly dependent on the interaction 
strengths. The new feature of discord offers promising appUcations in modern quantum technologies. 

A quantum many-body Hamiltonian H often possesses invariance under symmetry operations, exemplified 
by the particle number N and the angular momentum p, J^. Proper treatments of symmetry are of great 
importance in developing approximation methods for the Hamiltonian, in particular for strongly corre- 
lated systems. Unfortunately, ubiquitous approximation methods such as the mean field approach usually require 
breaking the symmetry of the Hamiltonian. For instance, the renowned BCS theory', proposed for supercon- 
ductivity and later for nuclear system, employs simple product wave functions of independent quasi-particles to 
reveal the underlying physics, at a price of violating the particle-number conservation. The entire system under- 
goes a phase transition to a symmetry-violating superfluid phase. In condensed matter physics, the particle- 
number fluctuation of a bulk superconductor is weak and the symmetry breaking does not play a significant role. 
On the contrary, the violation manifests itself in finite-size ensembles such as superconductive grains'' and nuclear 
systems. For instance, the nuclear properties obtained in the BCS treatment are the average of the target nucleus 
and its adjacent nuclei. The BCS theory successfully captures the dominant quantum correlation, i.e. the pairing 
effect between single particle (electron or nucleon) states and their time-reversals, but thoroughly washes out 
other correlations such as the correlations between different modes. These correlations might not be important 
for bulk superconductors but apparently exist in finite-size systems, and can be restored by the symmetry- 
restoration methods that go beyond the mean field approaches. 

In quantum information theory, quantum and classical correlations can be distinguished specifically. 
Quantum correlations characterize quantumness, measured by concurrence and quantum discord, of a correlated 
system. It should be of great interest how quantumness varies with the process from violation to restoration of 
symmetry. Here we study a two-level BCS model. The conventional BCS theory of the model violates the particle- 
number conservation. We shall restore the conservation by the projected BCS (PBCS) method, which projects out 
states with fixed particle numbers from the BCS wave function. It is noticeable that the lost quantumness is found 
when the symmetry is restored. Quantumness and symmetry are destroyed and reestablished simultaneously. We 
study the finite particle number effects on both concurrence and quantum discord in the case of large degeneracy. 
To our surprise, while entanglement drops greatly with the particle number from the maximum at two end points 
and then remains at half value, quantum discord behaves as an extensive quantity for a target system with a fixed 
degeneracy. This implies that the paring interaction is effectively transparent to discord. Physically, discord 
displays quantumness where a pair of qubits moves independently as if other qubits do not exist, whereas the 
existence of other qubits severely affects the entanglement. According to their definitions, quantum discord is the 
quantumness associated with quantum measurement, e.g. von Neumann measurement, on a system, while 
entanglement is from the density matrix itself of the system. The difference is crystal clear. Quantum discord 
seems to be a more general and fundamental nonclassical correlation than entanglement as pointed out in recent 
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literatures'"''. As proved in Ref. 6, the set of states with zero-discord 
has volume zero in the whole Hilbert space, very different from the 
set of non-entangled states which have finite volumes. Since almost 
all quantum states may be useful resources for discord-dependent 
tasks due to their positive discord, our result that discord acts as an 
extensive quantity, can be important in assessing availability of a 
quantum state in various quantum information tasks. 

The BCS theory' was proposed to determine the fully paired 
ground state of superconductors and later was translated into the 
theories of nuclear structure^ since the pairing correlation was found 
in low-lying states of even-even nuclei. In the past decade some 
authors studied the entanglement of superconductors in the frame 
of BCS model considering that it provides a solution to a quantum 
many-body problem with an explicit wave function. Ref 8 used 
concurrence to quantify the entanglement of the BCS state of super- 
conducting compounds. Ref 9 analyzed the finite-size properties of 
the two-level BCS model and discussed the entanglement properties 
of the ground state via concurrence. Without loss of generality, we 
now focus on the BCS application in nuclear systems, however, the 
conclusions obtained can be applied to any finite systems directly. 

The ground-state BCS wave function of an even-even nucleus 
reads'" 



|BCS)= n (uk + Vkalal)\0) 



k>0 



(1) 



which violates particle-number conservation. Here aj^ and a| are 
fermionic creation operators of single particle state k and its time- 
reversal k, respectively. Uf^ and Vj- are real variational parameters 
satisfying the normalization condition + = 1. These parameters 
are determined by variation of the energy, with the restrictions that 
the expectation value of particle number operator N is equal to the 
practical particle number N, i.e., (BCS|NjBCS) =N. 

The pure pairing interaction is widely used in nuclear BCS theory, 
for it provides a simple and powerful description of paring correla- 
tions in nuclei. A generic HamUtonian with the pairing interaction 
reads 



ki^alalaj.:ai^, (2) 



where e^-'s are single-particle energies and Gkk'kk' are the strengths of 
paring force dependent on the energy levels. The BCS theory is con- 
venient and widely used. However it violates the particle number 
conservation and the physics obtained from BCS treatment is an aver- 
age of the target nucleus and its adjacent nuclei. This symmetry break- 
ing can be restored by the PBCS theory", as introduced in Method. 

We now construct three bilinear fermionic operators, s'.^' = a|c!^ , 

S'*' = a^Uk and Sg''' = - ^fljfli + a^aj. — , which play the same role 

as spin operators (half of Pauli operators) and generate an su{2) 

algebra. The operator So has eigenvalues + - depending on whether 

the pair state [kjc) is occupied or not. We denote the empty and 
occupied with |0) and two states supporting a qubit. The equi- 
valence between fermonic pairs and qubits is discussed in Refs. 12- 
14. Specifically, a fermonic pair in (kjcj corresponds to the kth qubit. 



coupled in pairs. The BCS and PBCS wave functions are given by 
Eqs. (1) and (32) in Methods. 

Eq. (1) can be rewritten such that the BCS wave function becomes 
a direct product state of qubit states of different modes. Both con- 
currence and quantum discord are therefore vanishing. It shows that 
while violating particle number conservation, BCS wave function 
annihilates quantumness of different modes. 

Now we calculate the reduced density matrix of two qubits A and B 
under the PBCS wave function'". The Hilbert space of the two qubits 
is spanned by the computational base, {|00), |0l), 1 10), 1 11)}, where in 
each basis, the first (second) digit indicates the state of qubit A(B). 
There are three different types of reduced density matrices, deter- 
mined by the way qubits A and B occupy the two energy levels, and all 
of them are a particular case of two-qubitXstates'^"^"'"'^"^'' with Pi4 = 
p4i = 0, which in general reads 



(3) 



We distinguish the three types by superscripts (1), (2) and (3). Type 
( 1 ) denotes the case that both qubits A and B are in the lower level; (2) 
A and B in the lower and the higher levels, respectively; (3) both A 
and B in the higher level. The non-zero matrix elements of p'" are: 



/ Pn 
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0 
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P22 


Pa 


0 
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P32 


P33 


0 
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0 
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) 
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(1) (1) (1) (1) 
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mm{;),n,-l} / Q, -2 
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( — max{2,^ — Q2 



ill -2 

i~2 



^2 
p-i 

fi2 

p-i 

^2 
p-i 



,2\ P- 



,2\ P- 



(4) 



,2\ P-' 



where i^j^ j stands for the number of combination, the common 
factor C = uf^' u^' . For the second type, the matrix elements are 

™.|£;lo..n V ^ )U)[p-r)W2 

,=max{0,;,-n,} Vp~'~1 / 



(2) (2) / (2) (2) 
P23 = P\2 = V P22 P33 ' 



(2) 
P44 



_ ^ -""l^l-"'} / ill - 1 \ YY ^2 - 1 



2\P-' 



Results 

Consider a simplified pairing HamUtonian (2) (similar to the nuclear 
model in Ref 15, 16): two degenerate single particle energy levels, set 
as ei =0 and £2 = 1, with the degeneracies Qi and Sl2> respectively. 
The strengths of paring force are Gj-j-j.'^- ~ 4G,i(r = 1,2) for both k and 
fc' in the same level i and Gj.j.j.'p ~ 4Gi2 = 4G21 for k and k' in different 
levels. We focus on the low-lying states where all particles are 



The reduced density matrix of third type is similar to first one. When 
Qi = ^2 = fi> the two types are connected by a simple relationship, 
p'^^\p) = p'"(2n —p), due to the particle-hole symmetry with 2f2 — p 
being the hole number. 

For the X-type density matrix, one can derive the expression of 
concurrence as" 



C(/)) =max{0, 27p22P33-2VPiiP44}- 



(6) 
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Figure 1 | PBCS calculations for (a) Concurrence and (b) Quantum discord versus particle-pair number p at the paring strength G = 0.6, single 
particle energies = 0 and 62 = 1. Both levels have the same degeneracy SI = 20. The results of BCS are also drawn as a reference. AH connection lines 
between points are just for guiding eyes and the same below. 



The total correlation of the system can be written as'" 

J (p ) = P 1 1 l0g2 P 1 1 + P44 log2 P44 + (P22 + P33 ) log2 ( P22 + P33 ) 

- (Pu +P22) l0g2 (Pll + P22) - (P33 +P44) log2 (P33 + P44) (7) 

- (Pu +P33) l0g2 (Pll +P33) - (P22 +P44) log2 {Pl2+Pu)- 

The classical correlation is defined, based on the von Neumann 
measurements {Bj^} (fc = 0, 1) on subsystem B, as 

C(P)= -(Pll+P22)l0g2 (P1I+P22) 

- (P33 +P44) log2 (P33 + P44) -min{Si,S2}, 
where we follow the discussions of Ref. 20 and here 

Pll „ ,„„ P33 



(8) 



I P22 , P44 

-^^^^°^^p2;tp4:-^^*^°^^p2;tp4: 



and 
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1-0 
2 



, 1-0 1+0 

^^°S2 2 



1 + 0 



0 = 



(Pll + P22 - P33 -P44)^ +4P22P33 



(9) 
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Two-level case. We now consider a pair-correlated system with two 
degenerate levels. We first compute concurrence and quantum 
discord based on PBCS for a simple case that Q.i = Q.2 — Gn = 
G12 = G21 = G22 = G. The concurrence and quantum discord of 
three types versus pairing number p are shown in Fig. 1 (a) and (b), 
respectively. At first glance, we notice that both curves of 
concurrence and quantum discord exhibit symmetry about p = SI, 
perfectly for the second type, and approximately for the first and 



third types, as well as the particle-hole symmetry between first case 
and third case as mentioned above. Thus, in the following we only 
need to concentrate our discussion on the region 0 < p <Si and one 
of the three types, e.g., type 2. The more meaningful and essential 
result displayed in Fig. 1 is that concurrence and discord exhibit 
totally different behaviors. From Fig. 1 (a) we can see that 
concurrence starts with zero at p = 0 corresponding to the BCS 
ground state, jumps to the maximal value at p = 1, then decreases 
by near half and keeps almost unchanged with particle number. 
Clearly, it shows that symmetry restoration results in the 
reestablishment of entanglement. 

On the other hand. Fig. 1 (b) shows that quantum discord also 
starts with zero, indicating that the BCS state is not only unentangled 
but also classic-only correlated. The different behavior of discord 
occurs after p = 1, it does not drop like concurrence but increases 
with pair numbers until its highest value as the total maximal occu- 
pation number of the two levels are half-filled. This fact implies 
that concurrence and discord are different aspects of quantumness, 
owing to their distinct nature. There have been some reports about 
physical phenomena where entanglement and discord behave differ- 
ently. A typical example is that discord can indicate the Kosterlitz- 
Thouless phase transition in the XXZ model, whereas concurrence 
cannot"'■^'•^^ 

In Fig. 2 we plot together the curves of concurrence and discord, as 
well as the total correlation X(p) and classical correlation C(p) (refer 
to Ref. 22 for detailed expressions) vs particle number for the second 
type with Sli = Slj = 20. Apart from the remarkable difference 
between concurrence and quantum discord, we can also see that 
27 > C for the whole range of particle number, oppositing to the early 
speculation that C>'D for any quantum state^^"^'. 

The HamUtonian (2) can be exactly diagonalized. The numerical 
calculations show that the differences between PBCS and the exact 
solution wiU diminish with degeneracy and tend to vanish as the 
degeneracy is large enough. For a small degeneracy, e.g., Sli = SI2 
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Figure 2 | Discord, concurrence, total correlation, and classical correlation versus pair number p at the paring strength G = 0.6, and £1 = 20. Single 
particle energies are ei = 0 and 62 = 1. 



= 6, the fully restored quantum correlations measured by concurrence 
and discord in the exact solution are slightly larger than those in PBCS. 
It shows that PBCS is a perfect approximation to the exact solution in 
the sense of restoring the quantum correlations. However, superior to 
PBCS method, the exact diagonalization can even treat more compli- 
cated cases. In particular, we set Gn = 0.7, G12 = 0.6, G22 = 0.5 for the 
fixed Qi = Q2 = 6, and ili = 7, ^2 = 5 for the fixed Gn = G12 = G22 
= 0.6, respectively. Fig. 3 shows concurrence and quantum discord, 
and compare with a typical case that Qi = Q2 = 6 and Gn = G12 = 
G22 = 0.6. The deviations from the typical case are slight. 

Comparing Fig. 2 and Fig. 3, we also notice that for larger degen- 
eracy, quantum discord increases approximately linearly with par- 
ticle (hole) pair, when p,q<^Q, indicating that it acts as an extensive 
quantity like internal energy. We calculate the cases when Qi 7^ ^2. 
and Gn 7^ G22 = G33 G44, and confirm that these interesting 
properties remain the same. From Fig. 2, we can see that the quantum 
mutual information, i.e., the total correlation displays a much better 
linearity than discord. For p £ 8 the variation of total correlation 
versusp is very close to a straight line with slope ~ 0.42. For p £ 4 the 
quantum discord also coincides quite well with the straight line. An 
extensive quantity normally is the sum of the properties of separate 
noninteracting subsystems that compose the entire system. This 
implies astonishingly that the pairing interaction is effectively trans- 
parent to discord. For instance, discord as a function of the pairing 
strength G: 27(G)^1. 006(0.4), 1.003(0.5), 0.998(0.7), 0.997(0.8), 

0. 996(0.9) and 0.995(1.0) in unit of 27(G = 0.6) = 0.03168 when fi 
= 120 and p = 4. The values of X'(G) merely change with G, in 
particular when G is strong. Physically, discord displays quantum- 
ness where a pair of qubits moves independently as if other qubits do 
not exist, whereas the existence of other qubits severely affects the 
values of entanglement. Based on their definitions, discord reflects 
the quantumness from quantum measurement on the interested 
system, meanwhile entanglement is from the density matrix directly 
of the system. In comparison with entanglement, discord quantifies 
more general and more fundamental quantum correlations' ''. In 
particular, Ref. 6 shows that the set of states with vanishing discord, 

1. e., the set of classical states, has volume zero in the whole Hilbert 
space. This is very different from the set of separable (non-entangled) 
states, which exhibits finite volume. This result provides another 



perspective for understanding the differences between entanglement 
and discord. Considering that the existence of discord determines the 
nontrivial properties of quantum states, it was shown that discord 
may be the resource responsible for the quantum speedup"* in com- 
putational models. This implies that almost all quantum states are 
useful resources due to their positive discord. The feature, discord as 
an extensive quantity, is expected to have new applications in 
quantum information practices. 

One-level case. In order to understand the above results deeply, we 
consider a one-level system with degeneracy Q.. This is the limit of the 
two-level cases when level difference vanishes. The density matrix 
now reads 



fq{q-l) 
0 



fi(n-i) 



V 



0 0 

pq pq 
pq pq 



(12) 



0 0 p{p-l)J 



where q = Q — p is the hole-pair number in the single level. The 
strength of paring force G does not appear in the density matrix, 
which may be the reason why in the two-level case the entanglement 
and correlations are hardly affected by the strength of paring force. It 
turns out that it is paring correlation itself, rather than paring 
interaction, that determines these novel properties of a pair- 
correlated many-body system. Here we emphasize that p is 
symmetric under the interchange of p and q, or interexchange of 
states |00) and jll), so are the entanglement, classical correlations 
and quantum discord. In one-level case, all correlations have explicit 
analytical expressions, the concurrence is 

^^^^^ Q(n-i) [pq-VMp-i^q-i)]. (13) 

and the total correlation reads 
p(p-l 



fi(n-i) 

2pq 



log^|(^+^(^ 



log. 



fi(n-i) 

ipq 



n(n-i) °^n(n-i 



2p 



q[q-l) 



'n(n- 

Q 



(14) 
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Figure 3 | Concurrence and discord versus particle-pair number p from the exact solutions of Hamiltonian (2). Single particle energies are ei =0 and 

£2 = 1. 



As for classical correlation, we compare 
Si 



fi(Q-i) "'n-i fi(fi-i) '''n-i 
_ — P^log, ^ - 4^io& ^ 

I 02 o 1 r^/ O 1 A o2 ( 



(15) 



and 



with 



n(Q-i) ^'n-i fi(n-i) ''"n-i 



1-0, 1-0 1+0, 1+0 
S2 = -— log,^ r^°s^^ 



0= 



n(fi-i) 



(16) 



(17) 



We find that min {Si, S2} = S2, thus the classical correlation is given 
by 



C{P)- 



1-0, 1-0 1 + 0 



-log 



+ 



log; 



1 + 0 



^lo&^-^log,^ 



(18) 



and thereby the quantum discord is 



n(n-i) ^2Q(fi-i) Q(n-i) ^^q^q.j) 



2pq 2pq p p 

1-0^^ 1-0 i+Oj^ 1+0 



(19) 



These expressions are exactly symmetric under the interchange of p 
and q, therefore all of them are even functions oi p — q and their 



extrema lie at p = q = Q/2 and at the ends. In what follows we will 
focus our discussions on particle pair p due to the symmetry. All 
results are the same once replacing p with q. Concurrence C is a 
monotone decreasing function of p for p < Q/2 and reaches its 
minimum at p = Q/2. However, for a fixed degeneracy Q, the 
values of concurrence hardly vary with p, as shown in the above 
equations. Concurrence has the limit, as Q^l,p^l,p<CQ 



1 

Q-r 



(20) 



Q 



which is the same as its minimum at p = —. This indicates that 

^ 2 

concurrence behaves as an intensive quantity. On the other hand, 

2 

concurrence vanishes at large Q because of the factor — -— in 

^ n(Q-i) 

Eq. (13). This result may be used to choose materials as the resource 
of quantumness according to their degeneracy and the number of 
particle (or hole) pairs. 

On contrary to concurrence, the total correlation, classical cor- 
relation and quantum discord behave completely differently. They 
are monotonically increasing functions of p in the intervals p e (0, 
Q/2) with p = Q/2 as the maximum. The maximal values are 



Q-2 

max T ~l~ 



Q 



2Q 



max -52j^_i ' 2(Q - l) ^"^^ - 2 



1 /Q-2, Q-2 3Q-2 
log. 



4 \Q-1 



Q-1 Q-1 



log; 



3Q-2 
Q-1 



(21) 



-1 (22) 



1 /'3Q-2. 



'log2 



Q-2 



4 VQ-1 '"^^3Q-2 Q-1 
When Q^l, their limits are given by 



2Q 2Q 

log2 



= Q-2 



+ 1 



(23) 
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:max^-log23-l-0.189, 



-log^S^O.Sll, 



(24) 
(25) 
(26) 



It means that no matter how large the degeneracy is, the above 
correlations do not vanish but remain a constant. 

Under the large degeneracy limits with small numbers of pairs, the 
total correlation, classical correlation and quantum discord are 



Ap)\ 



2p 



C{P)\ 



log^, 



2p P , P 



(28) 



(29) 



It is thus clear that the total correlation increases linearly with par- 
ticle pair, while the quantum discord increases almost linearly, as 
observed in the two-level case where the slope is ~ 0.42 when = 40. 
We therefore further confirm, by the analytical one-level case that the 
conclusions of discord being an extensive quantity in the two -level 
case remains valid and should imply profound universality for any 
many-body system consisting of pairwise correlated particles. These 
theoretical results are obviously testable experimentally. For example 
strong correlated many-body nuclear systems may allow to choose a 
series of even-even isotones, which usually have the same single 
particle levels for valence neutrons, to examine the dependance of 
quantum correlations on the number of paired neutrons. We can also 
select a group of nuclei which have the same number of valence 
protons (or neutrons) but different valence shells to check our results 
for quantum correlations. 

Discussion 

We have studied symmetry restoration and quantumness reestab- 
lishment of the BCS theory, as well as the relationship between paring 
correlation and quantumness in pair-correlated many-body systems. 
Restored entanglement may not matter for bulk superconductors 
because it keeps on a low level with particle-pair numbers. 
Quantum discord, on the other hand, is an extensive quantity and 
grows linearly with pair numbers. From the perspective of paring 
interaction, a pair of qubits moves independently, while surrounding 
qubits strongly affects entanglement of the target qubit pair. The 
underlying origin of these effects should be from that discord is a 
measurement- dependent quantity, while entanglement is only deter- 
mined by the density matrix of the target system. Besides, because of 
the volume-zero of the set of classic states, these states can be hardly 
found in Hilbert space. Considering that in many physical systems, 
quantum discord behaves qualitatively similar to entanglement, it is 
significant to find physical entities where discord substantially differs 
from entanglement. An example is that discord signals KT phase 
transition while entanglement cannot. Our results display a new 
aspect of the differences, with respect to particle number effect, 
which may play a unique role in exploring strongly correlated sys- 
tems like high temperature superconductors since the quantity is 
hardly dependent on interactions in strong interaction regime'". 
The new feature may hint profound physics behind quantum discord 
and may offer promising applications for, e.g., certain quantum 
computations. 

It is worth mentioning that quantum discord can also be discussed 
in a more general setup of positive operator valued measurements 



(POVM)™"''. We can expect that, by considering POVM instead of 
von Neumann measurement, the main results will remain essentially 
the same. For X-states as considered here, extremization via POVM 
can be reduced to orthogonal projectors for a number of states^'. 

Methods 

BCS and PRCS. By introducing the Lagrange multiplier X, called the chemical 
potential, we determine the BCS parameters by variation of the BCS expectation value 
o{H' = H-XN 



(BCS|H'|BCS)=^(ei--^)v^+ ^ Gtm'vM. 

k>Q k.k'>0 

+ Gan-UkVkUfVk, 
k.k'>0 



(30) 



(27) and yields two quadratic equations for u^. and vf., 



/H^k 



(31) 



where 6^ — + ^ {Gkk'kk' + Gj:^k'kk')^k' ~ ^ ^^'^^ 8^P parameters 

Afc — — Gj^kkT ^k' Vk' ■ The BCS method provides the convenience that we can treat 

fc'>0 

a nucleus as a system of quasi-particles independently moving in a mean field. 
However it violates the particle number consei^vation. 

This symmetry breaking can be restored by projection techniques, of which the 

1 f 

method of residues'^ employs a projector — - — -q) ^ , , dz to act on the BCS wave 
function Eq.{l) 



(32) 



where C — z- andp — N/2 is the number of particle pairs. The integrand in the above 
equation is a Laurent series in C- Making use of the fermionic commutation relations 
for the operators a^, a^, the matrix elements of an observable can be expressed by the 
residues 



5^17! ) — - . 



dz 



2?:/ Jz(/'-^) + i k^k, 
For example, the expectation value of energy is written aS; 



n 



(T^IHI*") 



,K)-'\j:'kviR\{k) 



(33) 



(34) 



+ ^ [Gkkkkvlvl,Rl{k,k') + Gf^i^--^MkVkUk'Vt'Rl(k,k')] 

Kk'>0 

This projection method is termed as the PBCS theory. 



Concurrence and discord. Consider a mixed state p of two qubits A and B. By 
introducing operator p— ((7y(x)(7y) p* [<Ty(S)<Jy)^^, where Oy is they component of 
Pauli matrix, one defines the concurrence of p, 



C(p) = 



(35) 



where X., (/ = 1, 2, 3, 4) are the eigenvalues of pp in descending order. When C(p) > 0, 
qubits A and B are entangled. 

Quantum discord, as another kind of correlation, may exist even without entan- 
glement\ The quantum correlation features itself with many aspects, e.g., in char- 
acterizing quantum phase transitions^^. For the state p, quantum discord reads: 

V{p)=I(p)-C{p), (36) 

where X(p) is the quantum analogue of classical mutual information, defined as 

I{p)=S{p^)+S{p'>)-S{p), (37) 

with p"^'^' — TrB{A){p) denoting the reduced density matrix of the partition A (B), S{p) 
the corresponding von Neumann entropy. 2{p) is interpreted as a measure of total 
correlations in the composite system A + B, while C{p) is a measure of classical 
correlations, defined as 



C{p)=S{p'')- snp S{p\{B,}) 

m 



(38) 



where {B^} {A: — 0, 1) stands for the von Neumann measurements on subsystem B. 
From the above definitions, we can see that quantum discord is the quantumness 
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coming from quantum measurement, while entanglement is from the wave function 
(more generally density matrix) itself. 
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